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The invariant mass spectra of <f> — > K + K~ are measured in 12 GeV p + A reactions in order to 
search for the in-medium modification of <f> mesons. The observed K + K~ spectra are well reproduced 
by the relativistic Breit-Wigner function with a combinatorial background shape in three fij regions 
between 1.0 and 3.5. The nuclear mass number dependence of the yields of the K + K~ decay channel 
is compared to the simultaneously measured e + e~ decay channel for carbon and copper targets. We 
parameterize the production yields as a (A) — <jo x A a and obtain a^_^ K + K -— au_„,+ e - to be 
0.14 ± 0.12. Limits are obtained for the partial decay widths of the <f> mesons in nuclear matter. 

PACS numbers: 25.40.Ve, 14.40.Cs, 21.65. +f, 25.75.Dw 
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The properties of hadrons, such as mass, decay width, 
and branching ratio, have been extensively studied and 
well established in the history of particle physics. Recent 
interests have been extended to understanding how these 
properties are modified in hot or dense matter. This issue 
is of fundamental importance since such modifications 
can be related to the basic nature of QCD, spontaneous 
chiral symmetry breaking, i.e., the mechanism that cre- 
ates most of the hadron masses. Inspired by many theo- 
retical works related to this subject, several experiments 
including ours have been carried out. 

Amongst the many types of hadrons, the </> meson, 
which is a vector meson (1 ) of an almost pure ss state, 
has very attractive features for use as a probe to detect 
the possible changes in its properties. Its natural width is 
narrow (T^ = 4.26 MeV/c 2 ) [l| without any nearby res- 
onance; therefore, we may be able to clearly detect the 
possible mass modification. Theoretically, various mod- 
els predict the in-medium mass modification of the ^me- 
son both in dense and hot matter 0, 0, 0, 0, 0, 0, 0, 0] • 
The predicted decrease in mass at normal nuclear density 
is up to 40 MeV/c 2 and the width broadening is up to 45 
MeV/c 2 . Moreover, since m$ is only 32 (24) MeV/c 2 
greater than 2m^-± (2m^-o), the partial decay width 
^tp—tKK i s sensitive even to a small change in the spectral 
function of the <f> meson and/or kaon. Several theoretical 
models point out the possible change in the branching 
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There are few experimental reports on the search for 
the in-medium modification of the <f> meson. With regard 
to the ratio of the partial decay width T±_^ K % /r^j+j- , 
the experiments NA49 and NA50 at the CERN-SPS re- 
ported c/)-meson yields in the K + K~ and channels, 
respectively. There are discrepancies by factors ranging 



from 2 to 4 between these measurements 11, 12, Il3j |. 
Recently, the CERES experiment at the CERN-SPS re- 
ported new results in the (j> — > e + e~ and (j> — » K + K~ 
measurements [ijj- The results in the K + K~ channels 
arc in agreement with those of NA49, and both decay 
channels are consistent with each other. 

The present experiment E325, performed using the 
KEK 12 GeV proton synchrotron, recently reported the 
spectral modification of the (p meson in nuclear mat- 
ter measured in the e + e~ decay channel for the first 
time |l5(. We also reported the modification of p and/or 



ui mesons in Refs. jl6l Il7| . In the present study, new 



results are reported on the shape analysis for the <f> — > 
K + K~~ invariant mass spectra, and the nuclear mass 
number dependence of 0-meson production is compared 
between the e + e~ and K + K~ decay channels to de- 
termine whether the ratio of the partial decay width 
T^k+k- /rV_> e +e- depends on the nuclear size. This 
work is an advanced study over Ref. [ijj], where we re- 
ported the production cross-sections and their rapidity 
and transverse momentum dependences of the u> and (f> 
mesons measured in the present experiment. 

The cross section for a nuclear target of mass num- 
ber A is parameterized as a (A) = o-q x A a . When 
the <f> meson and/or kaon is modified in a medium and 
r^_ > ^+ Ar - /r^ e + e - changes, the ratio of the </)-mcson 
yield R = N^ K + K - /N^ e + e - becomes dependent on 
the mass number since a larger number of 4> mesons are to 
be modified in a larger nucleus. Consequently, by using 
two different nuclear targets A\ and A2, the difference 
in the a parameter between <f> — > e + e~ (a^ e + e -) and 
4> — » K + K~ (a^K+K-) can be related to R as follows: 
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FIG. 1: Invariant mass distributions of K + K~ pairs. The 
closed circles represent the observed distributions, and the 
solid lines represent the fit results with the expected (f> — > 
K + K~ shape on the combinatorial background shown with 
dashed lines. 
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It is important that most of the experimental effects are 
canceled in obtaining a. 

The details of the experiment can be found else- 
where [l9j]. A primary proton beam of 12 GeV with a 
typical intensity of 10 9 protons per 1.8-s spill was fo- 
cused on carbon and copper targets, and the spectrom- 
eter detected e + e~ and K + K~ pairs simultaneously. In 
the present article, we used the e + e~ -triggered data col- 
lected in 2001 and 2002 and the If +if "-triggered data 
collected in 2001. In 2001, one carbon (0.092 g/cm 2 ) and 
two copper targets (0.073 g/cm 2 each) were used, while in 
2002, one carbon (0.184 g/cm 2 ) and four copper targets 
(0.073 g/cm 2 each) were used. The targets were aligned 
in a line so that the same beam normalization could be 
used for all the targets. 

The invariant mass distributions of the K + K~ pairs 
are shown in Fig. [TJ To examine the modification in the 
mass shape as a function of /?7 (= p/m), the data were 
divided into three (3j regions. Each mass spectrum was 
fitted with a combinatorial background shape 24[ and a 
resonance shape of <fr — > K + K~ in the mass range be- 
tween 2mx± and 1.15 GeV/c 2 . The fits were performed 
by using the maximum likelihood method [20j | ; the "like- 
lihood chi-square x 2 " is shown in Fig.[TJ The combinato- 



- 


4 


(a) 


i 




ti ■ e+e 






X -«* 


¥ 




IF* 


if 


, i , , m^ud ^ 








events 











i (b) 

Ti 1 ' 

h 


200 


K 












h,, , *i„ 



rapidity 




rapidity 



FIG. 2: Observed kinematical distributions of (j> — > e + e~ and 
<fi — > K + K~ . Distributions of (a) f5"f, (b) rapidity, (c) pr, and 
(d) pr and rapidity. The points represent <j> — > e + e~, and the 
lines represent (f> — > K + K~ . The histograms for <j> — > K + K~ 
are scaled by a factor of ~3. 



rial background shape was evaluated by the event-mixing 
method. For the 0-meson resonance shape, the relativis- 
tic Breit-Wigner (RBW) distribution with the natural 
mass (too) and decay width (ro) was used after taking 
the spectrometer performance into account. The RBW 
distribution is given by 
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where 
,2/4 



TO 



to^) 1/2 m 



r(m) = r (q/qof (m /m) 
■ rn? K± ) 1/2 , and q Q = (m§/4 ■ 
The spectrometer performance was examined by a de- 
tailed detector simulation using the GEANT4 toolkit • 
The kinematical distributions of the <fr mesons used in the 
simulation were obtained from the nuclear cascade code 
JAM [23I, which reproduced our measured distributions 
fairly well [l8| . From the data, we obtained the peak po- 
sition to = 1019.43 ± 0.21 (stat) ± 0.04 (syst) MeV/c 2 
and the mass resolution a = 1.91 ± 0.29 (stat) ± 0.23 
(syst) McV/c 2 for the copper targets, while the simula- 
tion gave to = 1019.49 MeV/c 2 and a = 2.24 MeV/c 2 . 
The agreements arc significantly remarkable, and we have 
used the simulated mass shape for the fits in Fig. [T] 

The observed K + K~ spectra are well reproduced by 
the fit in all the /?7 bins. Therefore, the changes in mass 
spectra are not statistically significant in the <j> — ► K + K~ 
channel. Indeed, we have observed a mass shape mod- 
ification in the <p — > e + e~ channel fl5| in the very low 
/?7 region (^7 < 1.25) for the copper target data. In 
Fig. [21 we compare the acceptances of the K + K~ and 
e + e~ channels. In the region f3j < 1.25, we have very 
limited statistics for <\> — > K + K~, and we cannot ob- 
tain a reasonable fit for the K + K~ data. Thus, it is 
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TABLE I: 0-meson yields and the a parameters for the e + e~ 
and K + K~ decay channels evaluated by the fits. The first 
errors are statistical and the second are systematic. 



Pi 




JV (Cu) 


a 


(f) — > e~^~ € 
0-1.25 


331±32±6 


703±50±26 


1.18±0.07±0.02 


1.25-1.75 


702±51±9 


899±63±4 


0.87±0.06±0.02 


1.75-2.5 


945±62±14 


1096±71±12 


0.81±0.06±0.01 


2.5-3.5 


579±52±4 


626±56±5 


0.77±0.08±0.01 


4> -> K+K~ 








1-1.7 


99±20±18 


285±29±27 


1.38±0.13±0.07 


1.7-2.2 


143±24±19 


279±32±30 


1.14±0.12±0.03 


2.2-3.5 


177±26±11 


269±33±21 


0.99±0.12±0.03 



impossible to compare the 



K + K shape directly 



with the modified shape observed in the e + e~ channel. 
It can be concluded that in the region /?7 > 1.25, both 
<f> K + K~ and <f> e + e~ do not show signs of shape 
modification. 

We obtained the 0-meson yields for three bins for 
K + K~ and four bins for e + e~ as functions of j3~f. For 
the K + K~ decay channel, the </>-meson yields were ob- 
tained by integrating from 2m K ± to 1.07 GeV/c 2 after 
subtracting the background described before. To obtain 
the ^>-meson yields for the e + e~ decay channel, we em- 
ployed the same procedure as used in Ref. [l5j|; in that 
study, the spectra were fitted with a simulated c/> reso- 
nance shape and a quadratic background curve. The </>- 
meson yields were then obtained by integrating the data 
in the mass range between 0.9 and 1.1 GeV/c 2 after sub- 
tractin g th e background. Thus, the excess yield observed 
in Ref. [15j must also be counted as 4> mesons. 

Table |T] summarizes the yields for both decay channels 
with the systematic errors containing the contributions 
from the uncertainty in the background estimation [25| . 
In Table HI the a parameters are corrected for the target 
dependence of the experimental efficiencies. The most 
significant correction was the difference in the geometri- 
cal acceptance for each target position. The effect was 
estimated by the JAM simulation, and a changed by a 
maximum of 0.04. The other corrections due to the trig- 
ger, tracking, and analysis were estimated to be small 
and were included in the systematic errors together with 
those from the background estimations. In the analysis 
thereafter, the errors in a are the quadratic sums of the 
statistical and systematic errors. 

The observed a parameters are plotted as functions of 
/?7 in Fig. (Ha) and as functions of the rapidity (y) and 
transverse momentum (pr) m (b) and (c). The e + e~ and 
K + K~ decay channels cannot be compared directly be- 
cause of the difference in the detector acceptance between 
the e + e~ and K + K~ decay channels, as shown in Fig. [21 
Therefore, we determined a^_» e + e - two-dimensionally on 
the y-pT plane under the assumption that a^ e + e - is lin- 
early dependent on y and px, and estimated the values 




Si fi 



4. (b) 



rapidity 



(C) 



0.4 
0.2 




(d) 



Pt 



FIG. 3: a parameters of <f> —* e + e~ (closed circles) and (f) — > 
K + K~ (open circles) as functions of (a) f3y, (b) rapidity, 
and (c) pr- The horizontal error bars are the RMS values 
of the horizontal bins. The hatched bands show corrected 
a^,_ e + e - to be compared with a^ K + K - (see text). In (d), 
the differences between ot^_ K + K - and tt$_> B + e - in the kaon 
acceptance window are shown. The average value and error 
are also plotted as a hatched band. 



corresponding to the kaon acceptance windows [261 ]. To 
determine the plane a,^ e + e - (y,pr), the e + e~ data were 
divided into 3x3 bins in the y-pr plane, and they were 
fitted with the function a(y,px) = o, x y + b x p T + c. We 
obtained a = -0.32 ± 0.11, b = 0.13 ± 0.17, c = 1.24 
± 0.15, and x 2 /dof = 4.2/6, indicating that the above 
assumption is statistically acceptable. 

In Fig. [3ja), we show the estimated values of 
a cj>^e+e- {Pi) m the kaon acceptance window as a 
hatched band; this can be compared to the measured 
value of a^K+K - ■ The difference Aa in the kaon accep- 
tance is plotted in Fig. Eld). We expect that Aa should 
be zero when r<p^ K + K - /T < p_ te + e - does not change in 
nuclear media. Although it is interesting to observe that 
Aa increases when fii decreases, Aa at the lowest /3j 
bin is only 0.33 ± 0.17; the average value is 0.14 ± 0.12. 
Therefore, a,p^ e + e - and a,p^ K + K - are statistically the 
same in the measured kinematic region. 

On the basis of these results, the possible modification 
of the decay widths is discussed below. In a medium, we 
expect the total and partial decay widths T to change to 
r* according to the relation 



r;/r^ = 1 + hot (p/po) , 



Tl^K+K- T A^K+K~ = 1 + k K {p/po), (3) 



t>— >e+e- , 



1 + fc e (p/po) 



where T° is the value in vacuum and po is the normal 
nuclear density 



27|. We expect fc tot ~ fc^ since the 



meson mainly decays into KK as long as such decays are 
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FIG. 4: The obtained constraints on the in-medium modifica- 
tion of the partial decay widths of (f> — > K + K~ and <f> — > e + e~ 
(see text). 



kincmatically allowed. 

To estimate the relation between Act, kx, and k e , 
we performed the Monte-Carlo calculation described in 
Ref. [l5| . where mesons were uniformly produced in a 
nucleus and decayed according to the values of kx and 
k e . In the plane of kx and k e , the expected values of 
Aa were obtained. The results are shown in Fig. [5] as 
a colored contour, where the measured Aa and its error 
provide the 90% confidence limit. 

Next, we have investigated how the K + K~ spectra 
could provide constraints for kx- Since we observed a 
significant excess on the low-mass side of the 0-mcson 
peak in the <p — > e + e~ channel [l5| . we consider that 
a similar excess can exist in the <j> — > K + K~ spectra. 
We thus reanalyzed the <fi K + K~ spectra with the 
same fitting procedure as described earlier, except that 
the mass range from 2m K ± to 1.01 GeV/c 2 (0.01 GeV/c 2 
below the cf> peak) was excluded from the fit. Due to 
the statistical limitation, we combined all the kinematical 
bins into one. 

This procedure gives the amount of excess N ex as a sur- 
plus over the ^-meson peak and the background. Since 
the K + K~ acceptance is a function of the invariant mass, 
we corrected N cx /N,p and obtained 0.044 ± 0.037 (stat) 
± 0.058 (syst) and 0.076 ± 0.025 (stat) ± 0.043 (syst) 
for carbon and copper targets, respectively. Although 
these surpluses arc close to zero, they statistically limit 
the number of modified cf> mesons in the K + K~ spectra. 
The ratio N ex /N$ can be considered as the ratio of the 
number of (f> mesons decayed inside the nucleus to that 
outside it, i.e., N^/N^ = NT/NT*. 



The relation between N™/N^ nt and kx was also ob- 
tained by the Monte-Carlo calculation. We assumed that 
the decays inside the half-density radius of the Woods- 
Saxon distribution contribute to A^ n . By using the cal- 
culated relation, kx was obtained as 1.4 ± 1.1 (stat) ± 
2.1 (syst); this was the average value for carbon and cop- 
per targets. This result is plotted in Fig. [4j together with 
the upper-limit value, 6.0, at the 90% confidence level. 

The two constraints shown in Fig. [4] arc the first exper- 
imental limits assigned to the in-mcdia broadening of the 
partial decay widths of the <fi meson. To obtain the limits, 
we have renormalized the probability distribution func- 
tions by eliminating an unphysical region corresponding 

to r*/r° < 0. 

To conclude, the experiment KEK-PS E325 measured 
the </>-meson production via the e + e~ and K + K~ de- 
cay channels in 12 GeV p + A reactions. The observed 
K + K~ spectra are well reproduced by the relativis- 
tic Breit-Wigner distribution with a combinatorial back- 
ground shape. The nuclear mass number dependences of 
(f) — > e + e~ and <fr — * K + K~ are consistent. We have ob- 
tained the limits on the in-mcdia decay widths for both 
the decay channels. 
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